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. We found that the two-dimensional Schrodinger equation for 3 

electrons in an homogeneous magnetic field (perpendicular to the 
plane) and a parabolic scalar confinement potential (frequency ujq) has 
exact analytical solutions in the limit, where the expectation value of 
the center of mass vector R is small compared with the average dis- 
tance between the electrons. These analytical solutions exist only for 
certain discrete values of the effective frequency Co = \Jw% + (^) 2 - 
Further, for finite external fields, the total angular momenta must be 
Ml = 3m with m = integer, and spins have to be parallel. The 
analytically solvable states are always cusp states, and take the com- 
5-h ' ponents of higher Landau levels into account. These special analytical 

solutions for 3 particles and the exact solutions for 2 particles [ p^ j can 
be written in an unified form. The first set of solutions reads 

$ = Ui<k( r i - r k) m Pn,m(\ri - r k \) exp^-± u) n>m J2i rfj 
where p n ,m{x) are certain finite polynomials and w n>m is the spectrum 
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of the fields. The pair angular momentum m has to be an odd integer 
and the integer n defines the number of terms in the polynomials. For 
infinite solvable fields u>i there is a second set of the form 

where A a is the antisymmetrizer and the pair angular momenta 
can all be different integers. In both cases the first factor is a short- 
hand with the convention r m = r \ m \ e ima _ These formulae, when ad 
hoc generalized to N coordinates, can be discussed as an ansatz for 
the wave function of the N-particle system. This ansatz fulfills the 
following demands: it is exact for two particles and for 3 particles in 
the limit of small R and for the solvable external fields, and it is an 
eigenfuncton of the total orbital angular momentum. The Laughlin 
functions are special cases of this ansatz for infinite solvable fields and 
equal pair- angular- momenta. 

PACS classification: 

8.30.Vw (Quantum dots etc.), 73.30.Hm (Quantum Hall Effect) 
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1 Introduction 



Most work on correlated electron systems in a magnetic field (and a parabolic 
confinement potential) has been done adopting the following methods: Fi- 
nite particle number (N < 10) and finite field systems are tackled either by 
numerical expansion of the wave functions in antisymmetrized products of 
one- particle functions p -|| or analytical ad hoc approaches in the high field 
limit, where only the lowest Landau level (LLL) contributes 0-0. Other 
main streams are to use the Chern-Simons transformation Jl0| |~2[ and hop- 
ing that the transformed wave function can be guessed or approximated more 
easily, or to use models for the electron- electron interaction [|ll|] |12| (All the 
above references are mostly reviews). In this paper we are trying another ap- 
proach: We are looking for analytical solutions for few electron systems and 
trying to generalize them ad hoc to N particle systems. We use the genuine 
Coulombic electron- electron interaction and do not restrict ourselves to the 
lowest Landau level (LLL). In a previous paper fl3| it has been shown that 
for N = 2 (with Coulomb interaction between the electrons) there is a 'spec- 
trum' of discrete external field values for which the Schrodinger equation can 
be solved exactly and analytically. As shown below, these solutions comprise 
the Laughlin states for iV = 2 as special cases. The questions to be addressed 
in this paper are the following: Do similar exact solutions also exist for three 
electrons? If so, does the corresponding field spectrum agree with the spec- 
trum for N = 21 Is there any connection between the discrete field spectrum 
for solvability and the discrete fields (for given particle density) observed in 
the Quantum Hall effect? Are the Laughlin states still among these special 
states as special cases? To answer one of the questions in advance: We did 
not find exact analytical solutions for N = 3. However, if we consider the 
center of mass coordinate R = | J2f=i r i versus inter-particle distance as 
an small parameter and expand the Hamiltonian in a multi-pole series, the 
three- electron- system can be decomposed into 3 pair problems which have 
similar analytical solutions as the two- electron- system. The center of mass 
vector vanishes exactly in the classical ground state. Therefore, one should 
expect that our expansion works well for weak external fields or for systems, 
where after a Chern-Simons transformation and a proper mean field approx- 
imation (for finite systems!) the effective field is weak. Surprisingly, also for 
3 particles in the small R limit, the Laughlin states, which belong to infinite 
fields, are among the analytically solvable solutions. 
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The plan of this paper is the following. Sect. 2 gives a survey on the 
results of the exact solutions of the electron pair problem. Because in this 
paper the three- electron- problem is traced back to three pair problems, 
this seems to be helpful. In Sect. 3 we define an orthogonal transformation 
for the three- electron- problem which contains the center of mass R as a 
parameter, and then we expand the transformed Hamiltonian into a multi- 
pole series in R. In Sect. 4 it is shown that, in zero order in R, the transformed 
Schrodinger equation can be solved exactly and analytically for a certain set 
of external fields and total angular momenta. In Sect. 5, the eigenfunctions 
for 2 and 3 particles are written in an unified form and ad hoc generalized to 
arbitrary particle number. This expression is compared with the Laughlin 
and Jain states. In Sect. 6 it is shown that the analytically solvable states 
are just those states where a cusp appears in the energy versus total angular 
momentum curve. The accuracy of the multi-pole expansion is tested in 
Sect. 7 by calculation of the energy eigenvalues in first order perturbation 
theory in the dipole and the quadrupole term of the Hamiltonian. 
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2 Exact solutions for two electrons in an ho- 
mogeneous magnetic field 



In this section we summarize the results of a previous paper [13[ on the two- 
electron- problem and add some important subsequent unpublished findings. 
In particular, we add the asymptotic solutions to our former pattern, which 
had not been given the due attention in [pL3 1 , and incorporate the electron- 



electron coupling constant (3 explicitly, in order to be able to investigate 
the behavior of the exact solutions in varying the coupling strength. Com- 
pleting and reviewing the two- electron- problem is important, because in 
the present work the three- electron- problem is traced back to three two- 
electron- problems. It has been shown in ||T3| that the Schrodinger equa- 



tion for two electrons in an homogeneous magnetic field plus an external 
parabolic scalar potential has exact analytical solutions for a certain infinite, 
but discrete set of field values (hereafter referred to as 'solvable fields') []. 
Except for the asymptotic case of infinite external fields, which is part of this 
pattern, there is a one- to- one correspondence between exact solutions and 
solvable fields. Such solutions exist for singlet and triplet states as well as 
ground and excited states. A further qualitative feature is that these solu- 
tions occur, whenever a correlated state (with electron- electron interaction 
included) is degenerate with an uncorrelated one (without electron- electron 
interaction) Moreover, for each total spin and orbital angular momen- 



tum quantum number as well as for a given degree of excitation (ground 
state, first excited state, etc.), there is an infinite series of solvable fields 
which converges to zero. This means that the solvable field values are dense 
at zero. 

Now we are going to describe the general analytical form of the exact solu- 
tions. After introducing relative and center of mass coordinate 

r = r 2 -ri ; R=^(ri + r 2 ) (1) 

the Hamiltonian (in atomic units h = m = e = 1) 

H = E{^(p, + ^ A ^)) 2+ ^ 2 } + ^T^ + Hsvin (2) 



: If we speak of 'fields' without specification to a special one, we mean the effective 
oscillator frequency \ + (^) 2 , which is the relevant parameter. 
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decouples exactly. 

H = 2H r + ^H R + H spm (3) 



We follow the notation of [13j as long as not explicitly mentioned, u, 
is the oscillator frequency of the parabolic external confinement potential 
and A(r) = |B x r the vector potential of the external magnetic field. The 
center of mass degree of freedom behaves like a quasi-particle in rescaled 
external fields and the quasi-particle of the relative coordinate is a particle 
in rescaled external fields plus a rescaled repulsive Coulomb field originating 
from the e- e- interaction. The Schrodinger equation of the first problem is 
trivial, the latter problem, which will be considered here, is described by the 
Hamiltonian (see eq.(5) in [13] ) 



1 



P + -A r 

c 



1 2 



where f\to r = \uj 0) A r = |A(r). In polar coordinates r = (r, a), the following 
ansatz for the eigenfunction 

p imoL ll ( r \ 

■ m = 0,±l,±2,... (5) 



is justified, where the Pauli principle demands that m is even or odd in the 
singlet and triplet state, respectively. The Schrodinger equation H r 0(r) = 
e r 4>(r) gives rise to the radial Schrodinger equation for u(r) 

1 d 2 , lf_ 2 1\ 1 , 1_ 2 2 , 



2 dr 2 + 2 r " VV 2 + 2^" + 27/ M(r) = ir U{r) (6) 



where u r = ~a) = ly^ 2 , + (^) 2 , e r = I 6 = 6r — \ muJ c-, & c — ~ and 

oo 

the solution is subject to the normalization condition / <ir|w(r)| 2 = 1. In 

o 

dimensionless variables r — > v^V" an d e r — > e r /o) r the radial Schrodinger 



equation reads 

1 d 2 1 / 2 1\ 1 1 — 2 /3 



m - - + -(vovr) H — = — }tt(r 



2 rf(v^r) 2 2V 4/(v^;r) 2 2 V ry 



f ) ti(r) (7) 



2 The index V and 'i?' refers to the relative and cm. coordinate systems, respectively 
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The exactly solvable eigenf unctions have the following form 
where p(r) is a finite polynomial 

(n-l) 

p(r) (VE>)" (9) 

with n terms. The soluble fields and the corresponding eigenvalues are de- 
termined by the two requirements 

a n = ; a n+ i = (10) 

which guarantee truncation of the power series and therefore normalizability 
of the eigenfunctions. It is clear from (|7|) that the two truncation conditions 
depend from the parameters only in the combination I- and Conse- 
quently, we have two equations ( |T0"D for effectively two parameters. Tech- 
nically, we first calculate the solvable fields 4- from (16) in |T3] (with (3 
included) and then we get the corresponding eigenvalues from (17) in ||13|| . 
which we rewrite here in the form 

zr- = \m\ + n (11) 

Note that (|TT|) does not contain f3 explicitly, but only through the fact that 
the solvable u r depend on beta. The calculation of the solvable fields is 
non-trivial (and not repeated here in detail), but when they are found, the 
calculation of the corresponding eigenvalues is simple through (O) . Observe 
that the eigenvalues of (0) without e-e-interaction read |I8[ 



\m\+2k + l (12) 



where k is the node number. Comparison of (O) and (O) shows the above 
mentioned degeneracies of the interacting system with the noninteracting 
one. With n = 2k + 1, the interaction- free solutions fit into a generalized 
pattern of all analytical solutions. It can also be shown that equation (16) 
with (17) in fL3| , which defines the solvable |-, has always the solution 



g- = whenever n is odd. As an example, see (20a) in |T3| for the case 
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n = 3. (These infinite field solutions are not included in Tables 1 and 2 in 
Hl3| .) It is also clear from physical reasoning that for infinite external fields 
the e-e-interaction has no influence on the eigensolutions, because the kinetic 
energy dominates. Therefore, solutions without e-e-interaction are exact for 
Cj r — > oo. Sometimes we shall call these solutions 'asymptotic' solutions. 
The completed pattern has the following over- all- structure: For n — 1, 2 
there is one solution (ground states), for n = 3, 4 we have two solutions (one 
ground and one excited state), etc. (see Figure 1). Generally, the soluble 
fields are the smaller the larger the corresponding n is. For n — > oo the 
corresponding soluble uj t converge to 0. 

Now we consider the eigenfunctions. The series a v is defined by a two 
step recursion relation which reads for the soluble states (insert (17) into (14) 
in 0) 

ao = normalization constant 
1 P 
(2\m\ + 1) y/u r 

i r 1 p 2 oi u\ 

«2 = -71 i S —. — : : r 2(12— 1 J > On 

4(|m| + l)l(2|m| + l)w r V ! \ 

^ j=z a„-i + 2 {v - n - 1) 0^-2) (13) 



v [y + 2\m\) I v^V 

It produces rather complicated expressions for larger v. The recursion can 
also be started at the other end with normalization constant. 

The eigenfunctions of asymptotic solutions fit also into the generalized scheme 
For I- = 0, the recursion relation (p~3|) provides only non-vanishing coeffi- 
cients with even index, meaning, that p(r) is a function of r 2 . If we insert 
n — 2k + 1, and v = 2p with p = 0,1, 2, ... into ( |i"3"D we obtain the recursion 
relation 

(p — k — 1 
p{p + \m 

which belongs to the Generalized Laguerre polynomials. 

Pn=2k+i,m{r) = L^\uo r r 2 ) (15) 
3 We use the definition in Abramowicz, Stegun; Handbook of Mathematical Functions 



a p = ~i — n — FT a P-l ( 14 ) 
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Figure 1: Reduced energies (energy over effective oscillator frequency) ver- 
sus squared coupling constant over effective oscillator frequency for relative 
angular momentum \m\ = 1. The crosses indicate solvable states. The lines 
connect states with the same node number k. The termination index n is 
the same for all horizontal rows of crosses with the same ordinate. 



Consequently, the Generalized Laguerre polynomials are a special case of our 
polynomials p n ,m(j')- If we go in Fig.l along the line for ground states {k = 0) 
from the left to the right, then the polynomials for the exactly solvable cases 
have the following form. For n = 1 , p„=i,m.(r) = l)™\<jj r r 2 ) is a constant, 
for n = 2 (simplest case with finite fields) p n =2,m(j') is a linear function with- 
out a positive zero, for n = 3 we have a quadratic function without positive 
zeros, etc. Analogously, the exact solutions for the first excited state are all 
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polynomials with one node, but increasing order. 



As an overview, we give the formulae for the simplest exact solutions of 
the completed pattern (see also (19) and (20) in [|13j with j3 included). For 
n = 1 there is only the asymptotic solution 



p(r) 




1 



For n = 2 there is one finite- field solution 

(3 2 



p(r) 



2 (2 \m\ + 1) 

f3r 



1 + 



(2 \m\ + 1) 



Both former solutions are ground states. For n 
solution, which is a first excited state, 



= 



p(r) 



1 - 



cu r r 



(\m\ + 1) 

and one at finite fields, which is a ground state, 
P 2 



p(r) 



4 (4 \m\ +3) 

(3r 



1 + 



+ 



09 r) 



(2|m| + l) 2(2 \m\ + 1)(4 \m\ + 3) 



(16) 
(17) 

(18) 
(19) 

3 there is one asymptotic 

(20) 
(21) 



(22) 
(23) 



The corresponding energies follow from (O), and we put do = 1 for the con- 
stant term of the polynomial p(r) without loss of generality. 
In order to give an idea in what magnetic field range these exact solu- 
tions are located we consider the case without confinement (uo = 0). Then 
u r = uj c /4, where all solvable frequencies given above are in effective atomic 
units a.u.* (% = m* = (3 = 1). On the other hand, for GaAs we have 
B[Tesla] = ci; c [a.M.*]/0.1363. This means, that the largest finite solvable 
field ( for n = 2 and m = ) is u c = 2 a.u.* and B = 14.7 Tesla. 



10 



Now we add some words about the limit /3 — > 0. The wave functions of the 
asymptotic solutions (|- = 0) do not depend on (3 at all, indicating that 

they are robust against a variation of the e- e- interaction, f - = can 
be realized either by vanishing e- e- interaction (/3 — > 0) or infinite fields 
(uj r — > oo). On the other hand, the wave functions of the finite- field- so- 
lutions (It- = finite) evolve steadily into non- interacting ones for f3 — > 
indicating a strong dependence on the form of the e- e- interaction. 



If we insert for the center of mass system the ground state WF, then the 
total spatial WF in the particle coordinates reads (use (|5|), (||), and from fT3| 
formula (8)) 

$(r 1; r 2 ) = ( ri - r 2 ) m Pn,m{\^ - r 2 |) e"f ^> (24) 

where the first factor is a shorthand with the convention r m = r \ m \ e ma _ T n 
complex coordinates z = x + iy, z = x — iy, and with the opposite sign for 
the angular momentum [] fh = —in (compared with Laughlins notation [^] a 
'bar' has been added), this means 

r m = z m for m > (25) 
= z™ for m < (26) 



Therefore, for m < 0, ( p4[ ) reads in complex coordinates 

$(^, z 2 ) = (f! - ^ 2 )^ ^(Ifj - z 2 \) e -f (l*l a +W a ) (27) 

The solution for n = 1 (infinite field, p(x) = ao = const.) agrees exactly 
with the Laughlin- WF, in particular, fh = 1 is a determinant of two LLL 
functions, which corresponds to an uncorrelated full LLL. 
For m > we have 

z 2 ) = (^ - z 2 ) m VuAW - e-f (W a+ W a > (28) 



which is the complex conjugate of (p7|) and therefore the solution for the 
opposite direction of the magnetic field, if u = 0. In particular, n = 1 and 
m = 1 is a determinant of the two one- particle- states with k = and m = 



4 Observe, our m is the angular momentum quantum number itself and negative for the 
LLL. 
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and m — 1, in other words, it corresponds to the uncorrected solution of 
two full LLs. In summary, the solutions for n — 1 comprise the Laughlin 
WFs and the full LLs. It is tempting to assume that the other solutions 
(n = 2, 3, ...) have some relation to the other Quantum Hall states, but until 
now there is no prove for it. (see also Sect. 5) 
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3 Approximate decoupling for three electrons 

The Hamiltonian for three electrons in an homogeneous magnetic field (vector 
potential A(r)) and a harmonic scalar potential (oscillator frequency Uq) 
reads 

3 

1 = 1 l^K ' " ' 

3 

where H spin = g J2 s i ■ B. The goal of the following considerations is the 

i=l 

decoupling of the Hamiltonian into a sum of independent Hamiltonians for 
quasi- particles. 

For two electrons this happens automatically by introducing the center of 
mass and relative coordinate. This is mainly due to the peculiarity that 
there is only one interaction term r^f which contains only one of the new 
coordinates r and R. The peculiarity for three electrons, which can be taken 
advantage of, is the fact that the number of interaction terms is equal to the 
number of particles. This does not allow for an exact decoupling, but an 
approximate one into three noninteracting pairs plus a coupling term which 
is small in the strong correlation limit. 

It is easily shown that an orthogonal transformation f] leaves the kinetic en- 
ergy in a homogeneous magnetic field and the potential energy in an external 
harmonic scalar potential, invariant. On the other hand, the center of mass 
(cm.) of a classical system in the ground state vanishes and it is natural to 
assume that the cm. in the high correlation limit can be treated as a small 
expansion parameter. Therefore, we look for an orthogonal transformation 
which transforms the e- e- interaction in such a form that each term depends 
only on one of the new coordinates and the center of mass. Additionally, we 
demand that the center of mass is invariant under the transformation, which 
guarantees that if it is small in the original coordinates, so it is in the trans- 
formed ones. The transformation from the original coordinates r-j to the new 

5 Exactly speaking, for keeping the one- particle- terms of the Hamiltonian decoupled 
it suffices that the row vectors of the matrix in ( |30| ) are mutually orthogonal. However, 
using this more general type of transformation does not provide any advantage in our case, 
because the additional freedom destroys the symmetry among the particles. 
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2 T V - + c A < r - : 



i 2 2 

+ ^ 



+ E 



Th.\ 



+ H, 



spin 



(29) 



ones Xj, which fulfills all these requirements, reads 



" Xi " 




" 1/3 a b 




ri 


x 2 




b 1/3 a 






. X 3 . 




a b 1/3 






0=1/3- 


1/V3 and b = 1/3 


+ 1 


" ri " 




"1/3 6 a 




" Xi 


r 2 




a 1/3 6 




x 2 


. r 3 . 




6 a 1/3 




. X 3 



(30) 



1 / and its inverse is 



(31) 



x 3 



From ( |3~1"D it follows that 
ri - r 2 = v/3 (x - 
r 2 - r 3 = - Xi 

r 3 - r x = VS ( X - x 2 



(32) 



so that the Hamiltonian in the new coordinates is 

2 



H ■■ 
where X 



3 r 

E 

4=1 



T Vi + -A(xi 

l c 



1 2 2 1 







if 



spin 



(33) 



\/3 |xj - X| 

R is the center of mass in the new coordinates. It 



3 Zjj=l 

is possible (but complicated) to take care of the Pauli principle in the new 
coordinates. It is much easier first to transform the wave functions (WF) 
back to the original coordinates and then do the antisymmetrization . 
While being still exact, (|33|) is not exactly decoupled because X contains all 
coordinates. For X small compared with Xj , the e-e-interaction term can be 
expanded in a multi-pole series 



1=0 



where 



JL v 

V3 k 



1 



(34) 



(35) 
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In zero order in X, the Hamiltonian is decoupled and can be solved 
exactly and in many cases even analytically (see below). The general strategy 
for amending the zero order result is to consider the multi-pole terms in 
perturbation theory. As a response to a frequently asked questions, we want 
to emphasize the following. The approximation X = does not mean that 
the new coordinates Xj are no more independent. The one- particle part of 
the Hamiltonian is independent of X and still exact. X = only means 
making an approximation to the e- e- interaction term in H. 
It might be interesting to note that transformation fl30|), if applied to an 
Hamiltonian with an external Coulombic potential y (instead of the oscillator 
potential ^u^rf), transforms in zero order in X the e-e-interaction term and 
the potential energy in the external Coulombic potential into each other, i.e. 
leaves the Hamiltonian in lowest order virtually invariant. 



15 



4 Exact solution in zero order in the center 
of mass coordinate 



4.1 Pair Equation 

After expanding the kinetic energy, the Hamiltonian in zero order in X reads 
# (0) = I> + H spm (38) 

i=l 

with an effective pair Hamiltonian 

hi = + A + \wc k + -j= A (39) 

where = + (|cu c ) 2 , cj c = ^ is the cyclotron frequency, and U is the 
orbital angular momentum operator. This gives rise to the definition of an 
effective pair equation 

hi 0«(xi) = e q . g .(xi) (40) 

with the normalization condition 

d 2 Xi |0 ?i ( Xi )| 2 = l. (41) 



The subscript % at eigenvalues and eigenfunctions comprises all quantum 
numbers for the i th pair. In polar coordinates (x,a), the pair equation fl40" ) 
is satisfied by the ansatz 

; m = 0,±l,±2,... (42) 



must satisfy the radial pair equation 

^ d 2 1/ 2 1\ 1 1_ 2 2 ? 



, - m 2 - + -« ti(3;)=e!i(i) (43) 

2 dx 2 2 V 4/x 2 2 V3si V ; V ; V ; 

where e = e— \mio c and the solution is subject to the normalization condition 

oo 

/ dx \u{x)\ 2 = 1. In analogy to the two- electron- problem, we now use for 

o 

the radial part of the pair functions fl42|) the ansatz 

u( x ) = zH+g p(x) e -3 0x2 (44) 
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where p(x) is a polynomial, which is finite for solvable states. In this way we 
obtain for the pair function ( fi2"D in polar coordinates (x, a) 



(x) 



^7T 



(45) 



where x m is a shorthand for x ^e ima (see also Sect. 2). 



Because of the decoupling in zero order, the total eigenvalues and eigen- 
functions of can be obtained from 



E, 



91,92,93 c 9l ~ c 92 ' c 93 



spin 



91,92,93 l X l 5 X 2) X 3j 



'91 



v Xxj ■ (p q2 {X 2 ) ■ 



'93 



x 3 ) 



(46) 
(47) 



After inserting the transformation ( p2|) into the spatial part of the total 
WF (f47|), we obtain in the original coordinates 



91,92,93 



1 

i 



'9i. \ R- - "^(^ - >":•, 



'9-J \ ~^\ V? ' r ' 



(48) 



It is obvious that for any total WF the fields Cu q . = u q of the three pairs have 



to agree. If we insert (|4"5|) into ([4"8"D, we end up with 



91,92,93 



mi 



1 \ m2 

R- 7y|( r 3-ri)J p 92 

1 \ "13 

R - 7y^( r l - r 2)J P93 



R-^(r 3 -r i; 

R-^( ri -r 2 ; 



1 

2 ^9. 



0X1 .M -- ^IZ r i 



(49) 



In order to simplify the exponential factor to the present form we used the 
special property of our orthogonal transformation YU x i = Yi r f- This is still 
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the most general result. It holds for ground and excited state and can be 
simplified in special cases. For analytically solvable finite- field solutions the 
quantum numbers of all three pairs have to agree, because the solvable fields 
uj q for each pair depend on all quantum numbers contained in q. This implies 
in particular, that all angular momenta m ; in ((49|) agree. For asymptotic so- 
lutions, the solvable to does not depend on the angular momentum m and 
the node number k (or termination index n = Ik + 1) of the corresponding 
pair. Therefore, we can construct analytical solutions from pair states with 
different quantum numbers, i.e. different m« and ki. For asymptodic ground 
state solutions, which belong to rii = 1, fl49| ) can be simplified by omitting 
the polynomials, because they are constants. 

It should be mentioned that the eigenfunctions ( f47|) of form a com- 
plete set and therefore they can be used as a basis for a numerical solution 
of the full Schrodinger equation. Their advantage compared with basis func- 
tions constructed from one- particle states fU is that they contain 
an e- e- correlation cusp. As well known in Quantum Chemistry, the lack 
of this cusp in the basis functions gives rise to poor convergence in CI ex- 
pansions. In the present paper, however, we consider perturbation theory for 
improving the zero order result beyond the strong correlation limit R = 0. 



4.2 Pauli Principle 

The WF fl49|) does not yet fulfill the antisymmetry requirement. Taking care 
of this is particularly important for finding the exclusion principles for our 
type of WF. The question is: for what combination of quantum numbers 
and parameter values does the antisymmetrized WF vanish by permutation 
symmetry? In order to keep formula length under control, we introduce 
some shorthands and conventions. Firstly, the exponential factor in (^Uj) is 
fully symmetric with respect to permutations and therefore it can be simply 
omitted in the antisymmetrization procedure. Further, for the polynomial 
prefactor in ( fffi ) we introduce a single symbol 

*(x) = x m p(x) (50) 

Then the spatial total WF can be written as 

^W^i^rs) = T qitq2iq3 (r 1 ,r 2 ,r 3 )-exp(--Ug^rfj 
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T, 



91,92,93 



(ri,r 2 ,r 3 ) = t, 




(51) 



Now, only the polynomial prefactor T(r 1; r 2 , r 3 ) has to be antisymmetrized. 

The implementation of the antisymmetrization procedure and the an- 
tisymmetrized WFs can be found in the Appendix. The most important 
qualitative result is that a simple exclusion principle exists only for R = 0. 
It states that the WF vanishes by permutation symmetry: 
in quartet states S — §: 

if at least two pair functions agree and if the total orbital angular momentum 
M L is even. 

in doublet states S — \: 

if all three pair functions agree. 

Because for analytical finite- field- solutions all quantum numbers of the 
3 pairs have to agree, it follows from the Pauli principle (see above), that 
only quartet states (with parallel spins) can be given analytically. If m is 
the angular momentum of the pair functions in (|49| ) then the total orbital 
angular momentum in the solvable states is = 3m. In zero order on R, 
Pauli principle demands that m = odd. For finite R there is no restriction 
from the Pauli principle and m = integer. As will be shown in Section 6, 
it turns out that these solvable states are states where kinks in the curve E 
versus Ml occur and which bear the so called magic angular momenta (see 
Fig. 2). In short, only particular states can be solved analytically, but these 
states are the most interesting ones. 

For asymptodic solutions (with (3 2 /Co = 0) all quantum numbers can be dif- 
ferent, although the pair functions belong to the same external fields. This 
means that we can form total WFs from different pair functions and Pauli 
principle applies as given above. 
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Figure 2: The second term of the total energy in second order Taylor 
approximation (|76|) divided by c<i M\ /3 as a function of for the ground 
states (crosses). Crossed circles denote states which are forbidden by Pauli 
principle for R = (Laughlin approximation), but allowed for finite R. The 
states denoted by open circles are the ground states for R = if the states 
underneath are forbidden. The region Ml < 3 is omitted because it lies 
partly outside the scale and because it needs special consideration. 



4.3 Exact analytical solutions of the pair equation 

If we introduce in the radial equation for the relative coordinate in the two- 
electron- case (0) the same parameters as used for 3 particles (i.e. u and e), 
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we obtain 



Id 2 1/ j 1\1 l/c^= 2 V 2 (3} (N _ 2U , /e<*= 2 >\ v 



2 c/r 2 2 V 4/ r 2 2 V 2 J 2 r 
For avoiding confusion, the parameters and the WF in the two- particle- 
problem have been given the extra superscript N = 2. Comparison with 
the radial pair equation ( fI3| ) shows us that we obtain solutions of from 
solutions of (0) by simple rescaling. 
The 'spectrum' of soluble Co follows from 

u = I ^ (53) 
and the corresponding eigenvalues and eigenfunctions from 



2 



- e 
3 



,(N=2) (54) 



«*H5s u ™( T= ?s x ) (55) 

The prefactor in (|55|) has been chosen to retain normalization of the radial 
pair function, if u^ N=2 \r) is normalized. In any solution, the quotient 4 is 
equal in both problems and given by (see | fL3|1 ) 

~ ~(N=2) 

^ = w ^ ) =\m\+n (56) 

where m is the angular momentum and {n — 1) the highest power in the 
polynomial p(x). 

In this way we obtain from ([16]-|23]) and (|53|-|55|) the simplest exact solutions 
of (fl3"D as follows: 

For n = 1 there is only the asymptotic solution 

2 

= (57) 

10 

p(x) = 1 (58) 
For n = 2 there is one finite- field- solution: 

t = |(2H + 1) (59) 

" {x> = 1+ pf[Ti) (60) 
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Both former solutions are ground states. For n = 3 there is one asymptotic 
solution, which is a first excited state, 

2 

= (61) 

UJ 

= (hTT) (62) 

and one at finite fields, which is a ground state, 
2 

= 3(4|m|+3) (63) 



UJ 



4= Bx \^kP x 

v ( x ) = i + )H— + \ll L (64) 

P{ ' ' (2|m| + l) 2(2 |m| + 1)(4 |m| +3) 



The corresponding energies follow from (|5^D, and we put the normalization 
constant ao = 1 without loss of generality. The pattern of solvable states 
agrees qualitatively with Fig.l. Only the abscissa- values of solvable states 
are shifted. Those, who are not yet convinced in the correctness of these 
solutions are recommended to check them by insertion into (|43"1). 
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5 Comparison of our analytical solutions with 
Quantum Hall States 



In this section we confine ourselves to ground states and consider two cases 
separately. This will provide two different generalizations of the Laughlin 
WFs. We want to emphasize, however, that it would be possible (but not 
convenient) to include both cases in one formula. Secondly, we start with 
considering the case R = (zero order result) and add some words on the 
general case afterwards. Generally, the generalized formulae in this section 
comprise our analytical results for 2 particles and 3 particles. However, they 
could ad hoc be applied to any particle number and considered as trial func- 
tions. 



The first case comprises all solutions for three electrons with equal pair- 
angular- momenta in (|49| ) and the two- electron result (2^). In other words, 
it includes all finite- field solutions and and the asymptodic solutions with 
equal pair- angular- momenta. It can be written in the following unified 
form 

$ = II ( r * _ r k) m Pn,m{\*i ~ rjfc|) exj>(-- u> n , m ^2rf\ (65) 
i<k ^ 1 i ' 

It should be remembered that, apart from the case n — 1, the soluble field 
values u) n>m and the polynomials p n>m (x) depend on the particle number. 
Using complex coordinates as defined in Section 2, (|65| ) can be reformulated 

as 

$ = Yi(Zi - Z k ) m p n , m (\Zi - Z k \) Xl(^n,m) for m < 
i<k 

= complex conj. for m > (66) 

where Xi(^n,m) is a Slater determinant of LLL states for the effective fre- 
quency uj n) m and m — \m\ — 1. For m < our WF can also be rewritten 



using the Laughlin WF 0] 



$ = UPn,m(\Zi ~ Zk\) KTl «m) (67) 
i<k 

where the fields cD n>m occur in the exponential factors of the Laughlin function 
instead of the infinite field used in the original Laughlin function. 
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Our solutions ( |65| ) and the equivalent forms have the following properties: $ 
fulfills the Pauli principle, if m = odd (and m = even). It is an eigenfunction 
of the total angular momentum operator with eigenvalue Ml = m N( * N 2 ~ lS> , 
where N is the electron number. Apart from the asymptodic case n — 1, it 
has components in higher LLs due to the p(x)-factors. In the case n = 1 
(where p(x) = const and 4^ — 0) and m < 0, our WFs agree with the 
Laughlin states. 

<:r^i) = n(^r ex pf-^iEi^i 2 ) = n&-**)* xic^oces) 

For comparison, we also quote the Jain ansatz || 

= Plll J[(z t - z k r Xwifji) (69) 

i<k 

where - = ±— + m, Vlll is a projection operator onto the LLL and the 
determinant Xu*{&i) is for v* full LLs and taken at the asymptotically infinite 
frequency uj\. Apart from the special cases discussed above, (|66D and (|69|) 
do not seem to be fully equivalent. At least, both treatments contain the 
Laughlin WF as a special case, and there is the vague similarity that both 
have to do with higher LL components. 

The next property will be discussed using the formulae for iV = 3, but it 
would also apply to the corresponding generalized trial functions. If we go 
beyond the R = approximation, i.e. if we calculate the form of the pair 
functions in the transformed space Xj in zero order in X = R, but use the 
full back- transformation to the original coordinates Tj involving a finite R, 
then the solvable eigenfunctions are those in ( ^9|) with rri\ = m 2 = m 3 = m. 
Firstly, we have to remind that this function has to be antisymmetrized as 
discussed in the Appendix, because it is not automatically antisymmetric for 
m = odd as in the case R = 0. This provides a complicated expression. 
Due to this antisymmetrization it has simple zeros wherever two coordinates 
agree. However, there are no m- fold zeros as in the Laughlin WF, because 

the factors ^R — ^75 (r« — J do not vanish if two coordinates agree. This 

holds for all solutions (finite and infinite field solutions). This is a feature 
which agrees qualitatively with the Jain functions. 



24 



The second case differs from the first case only for N > 2. For the asymp- 
todic solutions in (|49|), different rrii are allowed, and for the ground state the 
polynomials can be omitted because they are constants. After introducing 
a new indexing for the angular momenta, which is more appropriate for the 



The second equation holds for < 0. For general m^, this function has to 
be antisymmetrized. For 3 particles, the antisymmetrized result is discussed 
in Section 4.2 and given in the Appendix. In this special case Pauli principle 
demands that = odd. $ in ([F(J) is an eigenfunction of the total orbital 
angular momentum with eigenvalue M L = J2i<k m ik, and lies completely 
within the LLL. It is apparent that the special case of equal = m = odd 
agrees with the Laughlin states. The changes produced by the extension of 
the considerations to finite R are analogous to the first case. 

A more detailed investigation of the applicability of the two generaliza- 
tions ( |65D and ([7(]) to iV-particle systems and in particular Quantum Hall 
states will be the aim of a separate work. 



iV-particle system, we obtain from ([49|) 




(70) 



(71) 
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6 Approximate analytical solution of the pair 
equation and magic angular momenta for 
three electrons 



In order to avoid misunderstandings, it should be told at the very beginning, 
why we are looking for approximate solutions, if there are exact ones. The 
answer is that the exact solutions exist only for certain fields and states, 
whereas the solutions of this subsection are for all parameters and all states. 
Although the exactly soluble states are the most interesting ones, the rest is 
necessary to prove certain cusp properties of the exactly solvable ones. From 
here on we consider only the case (3 = 1 and finite fields. 
As shown in []T3| | , sect. 3. 2, the effective potential in the radial pair equation 
d43f) can be expanded around its minimum into a Taylor series to second 
order. This cannot be accomplished fully analytically because the minimum 
position results from the zeros of a forth order polynomial equation. It is 
possible, however, to establish the approximate effective potential to order 
Tq 1 , where r^ 1 = (|u;) 1//3 is a small parameter in the strong correlation 
limit. The eigenvalues of the resulting oscillator equation can then be given 
analytically with the result 



e = c + ci m + c 2 m 2 + O(r 



where 



("a 

Cl 

c 2 



1 - 



UJr 



4 3 x /3 



3-1/3 -2/3 



where k is the node number. This provides a total energy of 
E = 3 c + ci M L + c 2 {m\ + m 2 . + ml) 



(72) 

(73) 
(74) 
(75) 

(76) 



where Ml = m\ + m 2 + 1TL3. It is clear that the Taylor expansion is the better 
the more symmetric the effective potential is. Therefore it gets poorer with 
increasing m and increasing Co. Nevertheless, this formula gives a qualitative 
understanding of the magic angular momenta. If we are interested in the 



26 



ground state for a given Ml, it is clear from (TTq) that it is formed by that 
set of rrii for which the sum of the squares of the rrii is minimal for a given 
sum of the This demand is met if the m; are 'as equal as possible'. As a 
example, for Ml = 2 there are two sets which provide the same Mi, namely 
(002) and (Oil), with the latter forming the ground state. It is also clear that 
the total momenta of the form Mi = 3m (multiple of three) play a special 
role because all three rrii can be equal in this case, but not for the other Mi. 
Fig. 1 shows the third (and only discontinuous) term of ( [76|) as a function 
of Mi. It is obvious that this curve has kinks whenever it is possible that 
all three rrii are equal, i.e. for Ml = 3m with m = 0, 1, 2, • • •. On the other 
hand, we learned, that equal rrii is a prerequisite for analytical solutions. 
Thus we conclude that the solvable states are always cusp states. 
Observe, that in the limit R = the Pauli principle demands that even m 
are forbidden (see Sect. 4.2). We also want to remind (see Sect. 5) that we 
obtain the Laughlin WF in the R = limit. If we go beyond this approxi- 
mation, the total WFs ( f48l) can be antisymmetrized for any values of the rrii. 
These facts elucidate the origin for the well known problem that exact diag- 
onalization procedures provide cusps also for those Mi, which correspond to 
even denominator Laughlin states [II]. This shows that the Laughlin ansatz (if 
applied to finite systems) has an additional symmetry (produced by putting 
R = 0), which is not present in exact solutions. These conclusions are not in 
contradiction with the symmetry considerations in [15 and |16|]. The latter 



papers find formula for the cusp (or magic) angular momenta, provided, the 
states under investigation are not forbidden by Pauli's principle. They do 
not have (and need) any explicit expression for the wave function. By the 
way, their general notion on the eigenstates is consistent with our small-R 
expansion. The harmonic approximation used in this section is also related 
to the harmonic approach used in []I7|, but not fully equivalent. However, 
our analytic approximations derived in Sect.s 3 and 4 is not harmonic. 
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7 First order perturbation theory and accu- 
racy of the strong correlation expansion 

Now we are going to calculate the contributions of the dipole and quadrupole 
term (|3^) and ( |37|) of the e-e-interaction to the total energy in first order | 
perturbation theory, i.e. 



(77) 



where the zero order wave function is generally given by (p7|) . For quartet 
states, where the total WF is just a product of spatial and spin part and 
for calculating matrix elements, the total WF in (|77"D can be replaced by 
the unsymmetrized spatial part. (It is simpler, however, to do the integra- 
tions in the transformed x^-coordinates rather than in the r,). The further 
calculation is straight forward and provides 



AE® 



1 



3\/3 
1 



k=l 



9^3 

where we defined moments 



k=l 

V - - T, E M i~ 3) Mi i ] 

' ' m k cy Z-^i Z-^i rn-h m k 

* k=ll(^k)=l 



Ml fc) 



dx x k [u m (x)Y 



(78) 
(79) 

(80) 



with u m (x) being the radial part defined in fl42| ) and given explicitly in special 
cases using ( |60|) and (0). If all three angular momenta agree: mi = mi = = m, 
the result simplifies to 



= — = 

v/3 m 

AE^ 1 



3^ 



jVf(-l) _ ]Vf(-3) M (2) 



(81) 
(82) 



It should be noted that the moment 3 ) , appearing in the second order 
contribution, diverges for m = 0. This is because for x — > the radial pair 



3 By first order we mean first order in all multi-pole corrections, but not first order in 



R 
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function goes as u m (x) — > x*-'™' 4 ^) and thus the integrand in (|8"0|) behaves 
like x~ 2 for small x. That is why the results for m = are missing in Table 2. 

For a test of the accuracy of our small-R expansion we use analytically 
solvable states only, i.e. quartet states with M L = 3m. We do this for 
magnetic field only (i.e. ujq = 0), because the confinement can be included 
afterwards by a simple rescaling of the parameters. For ujq = it follows 
from ([56]) and the definitions of Co and e that the zero order (in R) energy 
per electron reads 

E {0) = Q+ |m|) n 

3u c 2 2 1 j 

and, for comparison, the trivial result without e-e-interaction is 



T?(non-int) I , I |\ i 

= lm+l lj + fc + - (84) 

2>u c 2 2 V ; 

The formulae for the corrections in first and second order AE^ and AE^- 2 \ 



respectively, are given in (|81 ) and (0). Table 1 and 2 show the results for 
fixed m and varying n and fixed n and different m, respectively. E^ Taylor ^ is 
the result using the Taylor expansion of the effective potential in the radial 
Schrodinger equation as described in Sect. 6. Its agreement with the exact 
(analytical) solution of the radial Schrodinger equation E^ gives an account 
of the accuracy of the Taylor expansion. 

We conclude the following. 

i) Comparison of ^ non - mt ) an d E^ shows the contribution of the Coulomb 
interaction energy to the total energy. As to be expected, it grows with de- 
creasing u c and is tremendous for small u c (e.g. ten times larger than the 
kinetic energy for n — 10 i.e. u c ~ 10 -3 ). 

ii) For small \m\, the Taylor approximation provides a good tool for solv- 
ing the radial Schrodinger equation. For large \m\ the effective potential 
becomes so unsymmetric with respect to its minimum that its result goes 
fatally wrong. 

iii) The analytical zero order result E^°\ which is the main achievement of 
this paper, compares pretty well with the most precise result E^ 2 \ This is 
mainly due to the cancellation of AE^ and AE^ . The maximum error is 
of the order of 10 %. 
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Table 1: Total energy per electron in units of uo c for analytically solvable 
quartet states with total orbital angular momentum Ml = —3 [m = — 1) 
for u>o = (mag netic field only). E^ non - int ^ is the energy without electron- 
electron interaction, E <yTayloT ' > is the total energy where the effective potential 
in the radial pair equation is treated in second order Taylor expansion (ff6|), 
and is the exact result in zero order in R as given in fl8"3"|). AE^ 1 ' 2 ^ are 
the dipole and quadrupole contributions given in ([81]) and ([82]), and E^ is 
the sum of the three former contributions. The Zemann energy is omitted. 



n 




j^(non — int) 


^(Taylor) 
3u)r 


3uJr. 


aew 

3uJr. 


3uJr. 


3w r . 


2 


4_ 
9 


= 444.444 E-3 


0.5 


1.03796 


1 


0.154332 


-0.101713 


1.052619 


3 


2 

21 


=95.2381 E-3 


0.5 


1.49531 


1.5 


0.290373 


-0.126757 


1.663616 


5 


18.1896 E-3 


0.5 


2.49055 


2.5 








10 




2.20940 E-3 


0.5 


4.99381 


5 








15 




0.655360 E-3 


0.5 


7.49561 


7.5 









Table 2: Total energy per electron in units of uo c for analytically solvable 
quartet states for ujq = (magnetic field only) as a function of m with Ml = 
3m. (Only odd m are compatible with the Pauli principle.) All solutions 
belong to n — 2. The meaning of the other column heads is as in Table 1. 



m 




j^(non-int) 


j^iTaylor) 
3w r 


3uj c 


aew 

3lv c 


A£l 2 ) 
3uv 


EW 

3w r 


10 


4 

63 


=0.063 


10.5 


21.4662 


11 


0.164717 


-0.00834761 


11.15637 


3 


4 


=0.190 


3.5 


5.06036 


4 


0.161018 


-0.0283146 


4.132703 


2 


4 


=0.266 


2.5 


3.43552 


3 


0.158918 


-0.0435331 


3.115385 


1 


'I 


=0.444 


1.5 


2.03796 


2 


0.154332 


-0.101713 


2.052619 







= 1.333 


0.5 


0.892263 


1 


0.136400 






-1 




=0.444 


0.5 


1.03796 


1 


0.154332 


-0.101713 


1.052619 


-2 


4^ 


=0.266 


0.5 


1.43552 


1 


0.158918 


-0.0435331 


1.115385 


-3 


21 


=0.190 


0.5 


2.06036 


1 


0.161018 


-0.0283146 


1.132703 


-10 


4 

63 


=0.063 


0.5 


11.4662 


1 


0.164717 


-0.00834761 


1.156369 
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8 Summary and discussion 



We found that the two-dimensional Schrodinger equation for 3 electrons in 
an homogeneous magnetic field (perpendicular to the plane) and a parabolic 
scalar confinement potential has exact analytical solutions in the strong cor- 
relation limit, where the expectation value of the center of mass vector is 
small compared with the average distance between the electrons. These ana- 
lytical solutions exist only for certain discrete values of the external fields u). 
For finite external fields, analytical solvability demands that all three pair- 
angular- momenta agree what leads to total angular momenta Ml = 3m 
with m = integer. In zero order in R, Pauli principle allows equal pair- 
angular- momenta only for parallel spins and m = odd. The analytically 
solvable states are always cusp states, i.e. states where E(M L ) has a cusp. 
Further, these special analytical solutions for 3 particles and the exact an- 
alytical solutions for 2 particles could be written in a unified form, which 
contains only products over coordinate combinations ]7i<fc an d sums 
Conveniently, instead of using one formula we consider two cases fl55|) and 
(|70|). These formulae, when ad hoc generalized to N coordinates, can be 
discussed as ansatzes for the wave function of the N-particle system. These 
ansatzes fulfill the following demands: they are exact for two particles and 
for 3 particles in the limit of small R, and they are eigenfunction of the 
total orbital angular momentum. The Laughlin functions are special case, 
or in other words, both formulae provide two different generalizations of the 
Laughlin functions. 

Until now we know mainly that our WFs are analytically solvable states 
(within the approximations discussed above). It is also clear that states for 
an infinitesimally varied uj look quite different, i.e. the finite polynomials 
have to be replaced by polynomials with an infinite number of terms. For 
N=2 and 3 these 'neighboring' states are even explicitly known. It is not 
yet shown, however, if these special states have generally something to do 
with the Quantum Hall states (which show similar singular features), or if 
any physical quantity has any special feature in these soluble states. One 
encouraging fact is that the Laughlin states are special cases of our states. 
In prospect, it is possibly a good idea to look for similar exact analytical 
solutions in a spherical geometry instead of the disk geometry used here, 
because it is easier then to attribute a filling factor to each eigensolutions. If 
there is a connection between our states and Quantum Hall states, this would 
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imply some kind of inherent super-symmetry in the Quantum Hall states. 
Now we want to compare our treatment of three electrons with Laughlins 
0. Both approaches are approximate. Whereas he forms WF by antisym- 
metrization of one- particle states of the LLL, which is expected to be good 
for strong fields, we established an expansion, which is good in the strong 
correlation limit and which contains in general higher LL components. Con- 
sequently, we obtained a richer variety of solutions comprising the Laughlin 
states as special cases. 
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9 Appendix: Pauli Principle 



In order the obtain familiar looking formula, we rename T — > $ and t — > 0. 
The question here is, how the properly symmetrized spatial part $ has to be 
supplemented by an appropriate spin part in order to obtain a wave function 
which is eigenfunction of S 2 and S z (with quantum numbers S and Ms) and 
which satisfies the Pauli principle. S = Y.f=i s i is the total spin operator for 
all particles. For more than two particles this is a well established, but non- 
trivial procedure. The source of the difficulty is the fact that the spin space 
can be degenerate, i.e. there is more than one orthogonal spin eigenfunction 
Xi, (i = 1, .../) for given S and M (see Table 3 for N=3). 



Table 3: 

Standard spin eigenfunctions for N=3. a and (3 are the one-particle spin 
eigenfunctions for spin up and down, respectively. The first factor of a pro- 
duct of one-particle functions carries the spin variable T', the second factor 
carries '2' etc. 



s 


M s 


i 


Xi 


permutation sym. 


1 

2 


+1 

1 

2 


1 


[-{a[3+[3a)a+2 aa/3] 
^ [+{Pa+aP)P-2 P0a] 


symmetric for (12) 


1 

2 


1 

2 


2 


(a(3- (3a) a 
-j= 2 {(3a- a(3) (3 


antisymmetric for (12) 


3 
2 


+1 

l 

2 

3 
2 


1 


aaa 

[a{a(3 J r(3a) J r(3aa] 
^ \(3{(3a+a(3)+a(3(3] 
(3(3(3 


symmetric for all P 
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Generally, the total WF which fulfills our demands, can be calculated 
from (see e.g. [[19| ) 

2, 3) = Aa 2, 3) $(1, 2, 3) i = (1, ...,/) (85) 

with the antisymmetrizer 

A = 4ttE(-1) P ^ (86) 



VN\ v 

and P is a permutation operator. / is the dimension of the degenerate spin 
space (/ = 1 and 2 for S = § and |, respectively, for N=3) and (— l) p is 
the parity of the permutation V . The quantum numbers S and Ms of the 
spin part X^ as well as the quantum numbers of the spatial part $ are not 
indicated. The arguments (1, 2, 3) are spin or spatial coordinates depending 
on the function in question. Although being correct, (|8~5"D does not reveal the 
inherent permutation symmetry of the total WF. An equivalent symmetrized 
form is (see e.g. HQ 



/ 

m (87) 



= "75 E * * 

V / fc= i 

where the symmetrized spatial function is defined as 




** = \/7n£(-l) P tW)PS (88) 



and Uki(V) is an irreducible representation matrix of the permutation group 
for permutation "P given in Table 4. It is convenient to define column vectors 
$| (with i = 1,2) of the matrix Then the i th vector comprises all spatial 
information about the i th of the orthogonal states ^/j. 

Quartet S = | 

Because of / = 1 and because the spin eigenfunctions are symmetric against 
all permutations, the symmetrized spatial function is totally antisymmetric 
and we have 

* = X $ s ; <T = A$ (89) 

which is reminiscent of the permutation symmetry of the triplet state for 



N = 2. If we insert the solution ([51]) for $ we obtain a lengthy expression 
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Table 4: Irreducible representation matrices U(P) for N=3. For class and 
cyclic permutation symbols see [fl9 |. 



class 


P 


U(P) 


[13] 


£ 




1 
1 




[2,1] 


(1,2) 




1 

-1 




(2,3) 


l 

2 


-1 y/S 

Vs i 


(3,1) 


1 

2 


-1 -v^ 

-Vs i 


[3] 


(123) 


1 
2 


-1 x/3 
-v^ -1 


(132) 


1 

2 


-l -Vs 

V3 -1 
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which does not reveal anything. For R = 0, however, ([50]) implies m (— x) = 
(— l) m m (x) or in the short hand notation 

<l> m {i-k) = {-l) m <f> m (k-i) (90) 

with m the orbital angular momentum of the pair solution, we obtain 
1 

3 (3 - 1) 0i(l - 2) - (-1) M ^ 3 (1 " 2) 0i (3 - 2) 



(0i(2-3) 
+0 2 (2-3) 
+0s(2-3) 



2 (3 - 1) 3 (1 - 2) - (-l) A/i 2 (l - 2) 3 (3 - 1) 



0i (3 - 1) 0a(l - 2) - (-l) Mi 0i(l - 2) 2 (3 - 1) 



where Ml = Yh=i m i is the total orbital angular momentum. If we define a 
matrix of pair functions 



0i(2-3) 0i(3-l) 0i(l-2) 

02(2 - 3) 02(3- 1) 02(1 - 2) 

3 (2 - 3) 3 (3-l) 3 (l-2) 



(91) 



then our result can be written as a determinant or permanent[] of pair func- 
tions 



7S Det(s) 



even 



Perm(S) 



for M L 
ior Ml = odd 



(92) 
(93) 



From this fact we conclude the following important rule valid for R = 0: 

If at least two pair functions agree, Ml has to be odd. 

Otherwise the determinant vanishes. This means in particular, that the total 
ground state cannot be built up from identical ground state pair functions 
with m — 0, which would have the lowest energy without taking the Pauli 
principle into account. Therefore, the ground state consists of one pair with 
m = 1 and two pairs with m = and has a total angular momentum Ml = 1 
in the strong correlation limit. 

Observe the difference of these conclusions to the Pauli principle for one- 
electron orbitals. If $ were a product of functions which depend on one 



7 A permanent is similar to a determinant without the factors (— l) p 
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coordinate only (one-particle orbitals), which is equivalent to a Slater de- 
terminant for <£> s , in the quartet state all three spatial functions must be 
different, and the ground state would have a total angular momentum Ml = 
{m — 1) + m + (m + 1) = 3m where the integer m depends on the strength 
of the magnetic field. This holds if uj is so large that no excited one-particle 
orbitals are involved, i.e. in the weak correlation regime . Indeed, for strong 
fields the simple rule that Ml is a multiple of 3 is confirmed also by numerical 
calculations |ff|,IO|,|ffl. 



Doublet 3 = ± 

Because of / = 2 there are two degenerate and orthogonal functions which 
span a subspace. This level does not exist from symmetry reasons only if 
both functions vanish. We assume R = so that ( j90D holds and consider 
the two vectors one by one. 
For i = 1 and Ml = even we obtain 



fr 3 (l-2) D 12 (2-3,3-l) 
fr 3 (2-3) D 12 (3-l,l-2) 







1 









i " 




-1 


+ 2. 


-Vs 


1 




-1 


+ 2 


V3 



,(3-1) D 12 (l-2,2-3) 



where 



D 12 (2 



3,3 



Det 



<M2-3) 0i(3-l) 
2 (2-3) 2 (3-l) 



(94) 



and the other determinants are defined analogously: the subscripts refer to 
the quantum numbers of the pair functions involved and the arguments de- 
fine the arguments of the pair functions. 

For i = l and Ml = odd we obtain a similar formula as before but with the 

determinants in (^1) replaced by permanents called Pi 2 . 

For i = 2 and Ml = odd the result is also similar but with different column 

vectors 





' " 




1 



1-2) D 12 (2-3,3-l) 
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1 r y/3 



3) £> 12 (3-l,l-2) 



-1 



1 



-Vs 



0s(3 



1) £> 12 (l-2,2-3) 



} 



-1 



where D\ 2 is defined as in (|94|). 

For i = 2 and M L = even the determinants in the last formula have to be 
replaced by permanents. 

It is obvious that for <pi = <p 2 03 ) the determinants Di 2 vanish and 
consequently for Ml = even, vanishes and for Ml = odd, $2 vanishes, 
but both do not vanish simultaneously Therefore this state is allowed. How- 
ever, if all three pair functions agree, also those functions vanish in which 
permanents occur, because the sum of the prefactors sum up to zero. Con- 
sequently, the Pauli-principle tells us that for R = in the doublet state: 
All three pair functions must not agree. 

This is analogous to the one-particle model, where two orbitals may agree 
(if they occupy different spin states), but not all three of them. This means 
that, as for the quartet, the ground state consists of one pair with m = 1 
and two pairs with m = and it has a total angular momentum Ml = 1. 
Thus, in our limit and for vanishing magnetic field, the quartet and the dou- 
blet ground state energy agree. This degeneracy is lifted by a magnetic field 
because of the Zemann energy. 
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